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Within a quasiparticle framework, we reconsider the issue of computing the Bose-Einstein con- 
densation temperature (Tc) in a weakly non-ideal Bose gas. The main result of this and previous 
investigations is that Tc increases with the scattering length a, with the leading dependence being 
either linear or log-linear in a. The calculation of Tc reduces to that of computing the excitation 
spectrum near the transition. We report t wo ap proaches to regularizing the infrared divergence 
at the transition point. One leads to a a\/| In a|-like shift in Tc, and the other allows numerical 
calculations for the shift. 
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Inter-particle correlations in the weakly repulsive dilute Bose gas have been studied for more than half a century. 

' It is therefore quite remarkable that the dependence of the Bose-Einstein condensation temperature (Tc) on the 

, interaction strength has continued to generate controversy until this time. In this paper we present an analysis of the 

' longstanding Tc problem starting from a quasiparticle picture of the system in the un-condensed (high-temperature) 

^ ! side. 

• The model is described by the Hamiltonian 

a^: 

: ^ = E '^^i^^ +2v^ ^^^^ ^l-fk^l-k^^ , 

<N : , . . „^ . 

, where ek — k /2m is the free-gas spectrum, and b, b* are bosonic operators. In principle, one is interested in any 
possible interaction function f/(k), but the simplest possible form is a delta function in real space, so that C/(k) — U 
^ [ is momentum-independent. To first order U is related to the s-wave scattering length a hy U — Anh^a/m. A 
H ' dimensionless measure of the interaction is the quantity av}/^ , where n = N/V is the density. The program is to 
' ' examine the shift in the critical temperature Tc, as compared to a noninteracting gas of the same density, as a function 
y . of av}/^, for small an^/^. 

Q ' For several decades, the quasiparticle picture has been a central paradigm for studying correlated many-body 
^ , systems. Even in a highly correlated Bose liquid like liquid helium, a quasiparticle description accounts for the tran- 
^ sition temperature extremely well; e.g., the lambda point temperature can be obtained from the ideal gas expression 
•I— I , Tc'"' — (27r/m[C(|)]^/'^) rt^/"^ by replacing the helium atomic mass m by a quasiparticle effective mass, to*. The 
' effective mass in the liquid exceeds the bare mass, because of the inertia of the medium which has to make way for 
^ , any atom to move. Correspondingly, the transition temperature is lower for liquid *IIe than would be the case for an 
- - ■ ideal bose gas with the same bare mass. 

The question we would like to address is the following: how much information can one get about Tc from a 
quasiparticle description of the weakly interacting Bose gas? The simplest treatment, a Hartree correction, does not 
give any contribution to the transition temperature. The exchange contribution at mean field level (Fock), does not 
have any effect on Tc for a momentum-independent interaction, but gives a negative shift for a momentum-dependent 
interaction^". The Fock contribution to Tc depends on the details of the interaction, and for small interactions is 
weaker than the leading shift we are examining. The present understanding is that there is no mean field effect on 
the leading shift, and that higher-order correlations are required to calculate this leading shift in Tc. 

Interaction effects in the non-ideal Bose gas can be described in terms of the two-body t-matrix, i.e., the vacuum 
scattering amplitude. This quantity describes all possible collisions between two particles, without taking into account 
the fact that the surrounding medium has an effect on these collisions. Many-body corrections arising from the 
surrounding gas can be treated by using the many-body T-matrix, which is the sum of ladder diagrams. 

In this paper, we approach the calculation of Tc for a Bose gas using perturbation theory in the two-body t-matrix. 
This is to be regarded as a first step toward a calculation using the many-body T-matrix. The idea is that, if one 
can regard the system as a weakly interacting gas of quasiparticles, then perturbation theory in the quasiparticle 
interactions should give a reasonable starting point, and the effects of the surrounding medium can be treated as a 
correction to these results. 
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The work in this paper is not meant to be a definitive determination of for the weakly interacting bosc gas, but 
instead pursues the more hmited aim of calculating within a quasiparticle framework. As is typical for descriptions 
of critical-point properties, our approach faces infrared divergences. We elaborate on two ways to regularize these 
divergences. 

During the past several years, work by Baym, Laloe and collaborators^"®, has predicted a positive shift in Tc which, 
to leading order, is linear in the scattering length; ATc = aan^^^ + O{o?n^/^). Attempts to calculate the coefficient 
Ci, however, have continued to give fluctuating results. This has prompted the Baym group to predict^ a logarithmic 
contribution at the next order, i.e., a — a^lno contribution. The presence of such a non-analyticity might explain the 
difficulty in any numerical estimate of ci . 

Stoof, from a renormalization-group analysis^'^, has predicted an increase of the transition temperature proportional 
to a\ lna|, i.e., a non-analyticity at leading order. A leading non-analytic behavior would actually explain even better 
the widely varying results obtained in attempts to calculate the coefficient ci based on the assumption of linear shift 
ATc ~ cian^/^. One should also note that the work of Baym and collaborators relies heavily on power-counting and 
dimensional arguments. This kind of argument typically does not catch logarithmic corrections. The possibility of 
leading non-analytic dependence of ATc on the interaction, therefore, needs to be further examined. 

An older attack on the Tc problem is that of Kanno^^'^^. Using a "quasi-linear" canonical transformation, Kanno 
calculated the free energy and hence Tc of the Bose gas. Other derivations of the transition temperature include 
the canonical-ensemble calculation^^ of Wilkens et al leading to a negative shift, Schakel's effective-action theory^^'^® 
leading to a prediction for a positive linear shift, and Kleinert's recent five-loop calculation^'^. In addition, there have 
been spurious mean-field level predictions^*^'^'' for shifts in Tc proportional to the square root of the interaction. 

One could think of studying the Tc problem experimentally, e.g., in the context of trapped atomic Bose gases. 
Unfortunately, trap effects dominate in the dependence of Tc on the interaction for these systems. Size effects cause 
a decrease of Tc at mean-field level, far larger than the intrinsic effect under study in this paper. The "^He-Vycor 
system^^ might be more suitable for studying the Tc problem, but is yet to reach the accuracies necessary for identifying 
logarithm-like corrections. 



II. THE QUASIPARTICLE APPROACH 

The starting point of our analysis is the expansion of the single particle Green function to leading nontrivial order 
in quasiparticle interactions, namely: 



G(k,^) 



1 



^ - (ek - M + SR(k, Ei,)) - (I](k, z) - I]R(k, Ei,)) 
1 [E(k,^)-SR(k,£;k)] 



z-E], 



1 



du Ei(k, cj) 

TT {UJ — Z)(UJ — E\s) 



+... . (1) 



Here SR,i(k, lo) = d\t, rJmS(k, w-|-iO+) are the real and imaginary parts of the self energy, and i?k = ek — A' + S^nCk, £^k) 
is the quasiparticle energy. This expression has a number of attractive properties which would not appear in simple 
perturbation theory in the inter-particle interactions. 

In the quasi-particle approximation (1), the single-particle spectral function ^(k, w) = — 23mG(k, a; -h iO+) is given 
by: 



^(k,w) w 27r5(a;-^k) 



1 + 



duj' i;i(k,w') 



f dui' S] 
J ~W 



{oj' - EkY 
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Si(k,a;) 
{uj' - Ei^f • 



(2) 



Note that the form (2) automatically satisfies the spectral sum-rule, / du)A{\i,Lo) = 1. 

Wo proceed by writing the particle number iV as a momentum sum over the convolution of the spectral function 
and the Bose thermal distribution function, /(w) = l/{e^'^ — 1). 



N 



E 



div 
2^ 



^(k,w)/(a;) = No + V 



d^k 



/(i^k) + V 



d^k f dw[fiEk)^ f{u;)] 



r d-^k r du 
J (2^7 ^ 



(w - Ek^ 



Si(k,w) . (3) 



To proceed further we must resort to an explicit expression for the sclf-encrgy. We make the assumption that our 
quasi-particles interact weakly, so that we can work perturbatively. We limit ourselves to terms up to second order 
in the interaction. 
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FIG. 1: The two first-order and two second-order self-energy diagrams included in our calculation. Of the six topologically 
distinct second-order diagrams, only the two shown here axe necessary for momentum-independent interaction, because the 
other four can be obtained by Ist-order self-energy insertions in the first order diagrams. 



A. Perturbation Theory 

We first perform perturbation theory for the self-energy in the bare potential ?7(k), retaining explicit momentum- 
dependence. The diagrams that need to be calculated are shown in figure 1. After performing the Matsubara sums, 
we get: 



d'P f [/(gq) + /(gq)/(gp+q) + ./(gq)./(gp-k) ^ /(gp+q)/(gp-l.)] 



(27r)3 J (27r)3 z - (ep+q + ip-u - gq) ^ 

(7(p) [t/(-p) + C/(k - p - q)] . (4) 

Here Ck = Ck — A*- For a momentum-independent ?7(k) = U, the first-order Hartree and Fock terms give equal 

contributions, as do the two second-order diagrams. However, Using U(k) = U directly in (4) causes ultraviolet (UV) 
problems in the second-order self-energy, due to the term with a single Bosc function. The prescription for removal of 
such divergences is to shift to a description in terms of the vacuum scattering amplitude t, which is a more physical 
quantity than the bare potential. The bare potential C/(k) is related to the t-matrix by 

7T^^.^ ^. ^^.\t\ \ f '^"^^ [(p - k, q + k |f | p - k^ g + k') (p - k', q + k^ \i\ p, q)] 

C/(k) =aie(p-k,q-Fk t p,q) - / 7^-^ '—^ (5) 

J (/iTT)'^ Cp + Eq - Ep-k - Sq+k 

When eq (5) is used to replace the bare potential by the i- matrix in eq (4), we get an 0{t^) term from the first-order 
self-energy that exactly cancels the UV-divergent term appearing at second order. There is now no problem in using 
a low- momentum limit in which the i- matrix reduces to a momentum-independent constant t = Anh'^a/m. The real 
(on-shell) and imaginary parts of the self-energy now reduce to: 

EH(k,e-k) « 2tn + 2t^ f J% [ ^ fi^^^fi^^^^) + /(^~.)/(gp-^) - /(gp+.)/(gp-k) 

J (27r)3 J (27r)3 e^ + Ck- Cp+q - Cp-k 

and 

/d'^P f d^ Q 
(2^ J (2^ ^ ^'^ " + ^P"'' " ^''^^ 

X [/(gq) + /(gq)/(gp+q) + /(gq)/(gp-k) - /(gp+q)/(gp-k)] • (7) 



B. Self-consistent Perturbation Theory 



We will use (6) and (7) to define a self-consistent quasi-particle approximation, formally second order in the vacuum 
scattering amplitude. This is done by replacing the free-particle dispersion functions (cq's) on the right hand sides of 
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(6), (7) by the quasiparticle dispersions, £^q's. The imaginary part now becomes 

X [f{E^) + /(i^q)/(^^P+q) + /(i^q)/(J?p-k) - f {E^+^) f {E^_^)] . (8) 

When this expression for the imaginary part is used in eq (3), the Sj contribution cancels out exactly, resulting in 

i.e., the number of quasiparticlcs is equal to the number of particles. Since we will be using these results above and at 
the transition, we have neglected the zero-momentum occupancy Nq as compared to A'^. We note that the imaginary 
part of the self-energy, i.e., effects of quasiparticle broadening, has dropped out of the expression for particle number 
and hence does not play any role in our determination of the transition temperature. 

Finally, making the real part of the self energy (eq 6) self-consistent, we get the following expression for the 
quasi-particle energy: 



= Ck - /U 



+ i]R(k,i;k) 



d^P f d^q /(iJq)/(£;p+q) + /(iJq)/(£;p_k) - /(i5p+q)/(i^p-k) 



f d^P f 

J {2nrJ 



(27r)3 



-£^q + -E^k — -E'p-l-q — Ep-u. 



(10) 



III. CRITICAL TEMPERATURE AND SPECTRUM 



To derive an expression for the change in the critical temperature due to quasi-particle interactions, we first consider 
the density of an ideal Bose gas, both at its transition temperature Tc^^ = (27r/m[C(|)]2/^) n^/^, and at temperature 
Tc which is the transition point of the interacting Bose gas. 

Here Tc = T^°^ + AT^. We are looking for the difference of transition temperatures of an interacting and a non- 
interacting gas at the same density, therefore n^°\Tc^^) = n(Tc) = n. Evaluating (9) at Tc and subtracting equation 
(11) we obtain the leading correction to Tc due to interactions: 

ATc 2 /■ d^fc 

-3;^/(2;o^[^^(^^)-^^(^^)]' (12) 

Here fc{x) = l/(e^/'^-^= — 1) is the Bose distribution fimction evaluated at the critical temperature; and ^k = E-\^{Tc) = 
ek — jJ-iTc) + SR,(k, ^k) = Ek + SR,(k, ^k) — Sr(0,0) represents the single particle excitation spectrum at T^. The 
problem of calculating the relative shift has thus been reduced to calculating the spectrum at the critical point, ^k- 

At this stage, we are able to present a simple argument for the increase of Tc due to the addition of a repulsive 
interaction. This is important in light of the continuing appearance of negative-shift predictions in the literature^'^'^^. 
The major contribution to eq (12) comes from the infrared; therefore approximately 

/•cutoff 

ATc DC / dk k^[e^^ - ^-^] . 
Jo 

If the quasiparticle spectrum ^k is "harder" than the bare spectrum ek = k^/2m, i.e., if ^k is sub-quadratic in k or 
quadratic with an effective mass m* < m, then the ek integral dominates and the shift is positive. On the other hand 
if the spectrum at Tc were "softened" by the interaction, one would have a decrease in the transition temperature. 
Since the spectrum ^k at Tc is known to be intermediate between linear and quadratic, the weakly interacting Bose 
gas has a positive shift in the transition temperature with the introduction of a weak repulsive interaction. 
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A. Spectrum from Perturbation Theory 

It follows from (10) that satisfies the nonlinear integral equation, 

Ck = Ck + 



J {27TfJ (27r)3 ^p + Ck - Cq - Cp + k-q 

.,2 f d^P f /(^p)/(^q) + /(Cp)/(^p-q)-/(^q)/(^p-q) 

y (27r)3 7 (27r)3 _ _ i^^^ 

In the perturbative quasiparticle approach, equations (12) and (13) completely determine ATc- 

If the quasiparticle dispersion were quadratic at small momenta, the self-energy integral in eq (13) would be 
divergent in the infrared. A power-counting of this equation shows that, for self-consistency, the dispersion at small 
momenta should be sub-quadratic, ^k ^ fc^/^, modulo possible logarithmic factors which power-counting cannot 
predict. The power-counting arguments were first devised by Patashiskii and Pokrovskii^°. 

Unfortunately, this approach using second-order perturbation theory over-modifies the infrared spectrum. The 
long- wavelength spectrum at the transition point should be ^k ^ fc^^**, where rj is the anomalous dimension. The 
present system falls in the same universality class as the N = 2 quantum rotor model. For this universality class, a 
recent Monte Carlo calculation^^ has produced the value r] « 0.038. In contrast, eq (13) gives r] = 0.5. 

Acknowledging that the present framework over-modifies the spectrum at Tc, we proceed to analyze eqs (12) and 
(13) in order to extract the shift in the critical temperature due to interactions. 



B. Infrcired Divergence 

There are several ways of dealing with the infrared divergence. One is to treat cqiiation (13) sclf-consistently and 
determine the self-consistent ^(k) numerically; this procedure is outlined in section V. This approach has also been 
explored by Baym and collaborators. 

In the next section, we take a difi^erent approach, motivated by the work of Kanno^^. We make an expansion in 
5]R(k,^k) that allows us to use the chemical potential at the critical point, fi{Tc) = iic, as an infrared cutoff. This 
approach allows us to express the relative shift in the transition temperature as an expansion in /9cMc = Mc/^b7c, 
wh ich its elf can in turn be related to the scattering length. This leads to a transition temperature shift ATc that is 
a^JK.{a) to leading order, where K.{x) is a function approximately like | ln(a;)|, to be specified in the next section. The 
appropriateness of this procedure is difficult to gauge, but it provides strong support to the possibility of non-analytic 
corrections to the transition temperature at leading order. 

Another infrared regularization procedure that has appeared in the literature is A.M.J. Schakel's prescription of 
analytic continuation^^. The connection between Schakel's treatment and other calculations awaits further study. 



IV. NON-SELFCONSISTENT CALCULATION 



We expand /c(^k) = /c(ek - Mc + ^k{K Ck)) around ek - Mc: 

Um - /c(ek - Mc) + /c(ek - Mc)SR(k, ek - Mc)) • (14) 

This provides us with an infrared cutoff, — /Xc, for our integrals. The approximation is non-selfconsistent in the sense 
that a modification of the spectrum to a non-quadratic form no longer plays any role. Using this expansion, the 
relative shift in Tc becomes 

^ " 'Llwf ^^^^"^ " ''^^ ~ ^^^"'^^ " ^ / (^^^^'' " " ^^^^ 

Both integrals can be expanded in |/3c/Uc|, as detailed in appendix A, each integral being 0(-\/|/9cMc|) at leading order. 
The result is 

+ o{M = I [c(i)] VMM + o{M (16) 



ATc 
t(0) 
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3n 



A3 
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FIG. 2: Approximation for eq (17). Solid curve is exact numeric solution; dashed curve shows approximation l/^cA^cI ~ 
167r(a/Ac)^| ln(a/Ac)|. Inset shows comparison for the relative shift of Tc; solid curve is numeric solution and dashed curve is 
the ~ 0\/| ln(a)| approximation. 



The infrared eutofF needs to be expressed in terms of the scattering length; this is achieved by using the condition for 
the transition, /Xc = Sr(0, 0). The relation is 

/3cMc « ln(-/3c/"c) , (17) 

also derived in appendix A. The functional dependence of |/3c/ic| on the scattering length is thus similar to ~ 
(a/Ac)^| ln(a/Ac)|, and therefore we get a shift approximately of the form ATc ^ a-^/\ lna|. 
The approximation of eq (17) by |/?cMc| ~ 167r(a/Ac)^| ln(a/Ac)| is displayed in figure 2. 



V. SELF-CONSISTENT APPROACH: NUMERICS 



A self-consistent treatment of (13) involves a calculation of the (sub-quadratic) quasiparticle dispersion ^(k). The 
principal values in (13) make it difficult to treat numerically. A "high-temperature" approximation, obtained by using 
fiO ~ produces a simpler version without the zeroes of the principal value form: 

a = ek-2t^/^/^ . (18) 

While the same dimensional arguments apply to (18), it is likely to lose information about logarithm-like factors. 

Since the self-consistent spectrum has infrared behavior ^(k) ^ fc^^^, we look for solutions of the form ^(k) = 
k^/'^(j{k) + k"^ where a{x) is some smooth function that decays to zero for large x, and is constant or logarithm-like 
at small x. (Here the momenta and ^ have been rescaled to be dimensionless, by factors ^J(i/2m and /3 respectively). 

Figure 3 shows a numerically determined spectrum ^(k). Once ^(k) has been computed for a particular a/ A, the 
relative shift in Tc can be computed using our basic eq (12). The AT^ vs. a plot in the inset shows a definite linear-like 
increase, as opposed to, say, a \/a or a^f^ dependence. However the existence, or lack thereof, of a logarithm-like 
factor is difiicult to prove from numeric data. 



VI. CONCLUDING REMARKS 



To summarize, we have outlined the calculation of Tc from a perturbative expansion in the two-body f-matrix, or 
equivalently the scattering length a. The transition temperature is expressed in terms of the quasiparticle spectrum 
at Tc, and hence the problem is transformed to one of calculating the critical-point spectrum = Sk(Tc). 

In developing the formalism, we find that the use of the two-body f-matrix removes the ultraviolet divergences 
that appear in a perturbative calculation in terms of a momentum-independent bare potential U . We also find that 
the width of the quasiparticle peak (Sj) does not afii'ect Tc. Whether this is true at all orders, or a peculiarity of 
second-order, and whether this is true for a more accurate treatment in terms of the many-body T-matrix, remain 
open questions. 
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FIG. 3: Quasiparticlc spectrum, for (a/A) = 0.0005. (The "fit" indicates quality of iteration - line shows penultimate, and dots 
show final, iteration). Crossover from k'^^'^ to behavior is clear; dotted line is the noninteracting k^ spectrum for comparison. 
Upper inset shows the crossover function cr(fe) for several o/A values (9 x 10~®, 8 x 10"'' and 4 x 10~^). Lower inset plots 
relative shift of Tc against interaction. 

The new feature of this work is an unusual approach to deahng with the infrared divergence that appears in 
treating the spec trum at Tc. This regularization scheme bypasses the issue of spectrum modification, and results in 
a non-analytic (Vln a-like) correction factor to the linear shift of the transition temperature. 

The major remaining issue in the Tc problem seems to be the possible existence of a logarithm-like factor in the 
leading a-dependencc. The resolution of this problem should lie in a clear treatment of the spectrum at Tc. Of the two 
major approaches, the perturbative approach alters the quadratic spectrum too strongly, while the RG calculation^ 
neglects spectrum alteration. 

The reason for perturbation theory over-estimating the spectrum modification remains to be examined in detail. It 
has been suggested® that the reason is that perturbation in a or t does not take into account the fact that the effective 
interaction strength flows to zero under renormalization group fiow, at the critical point. An alternate way to view 
this is to note that the many-body T-matrix (sum of ladder diagrams) vanishcs''''''^^ at the transition temperature. 
A perturbative calculation in terms of the many-body T-matrix, as opposed to one in terms of the two-body f-matrix 
that we have considered here, may therefore be expected to give a more realistic modification of the quasiparticle 
spectrum at Tc. A calculation of in terms of the many-body T-matrix has not yet appeared in the literature. 
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APPENDIX A: DETAILS FOR THE NON-SELFCONSISTENT APPROACH 

In this appendix we evaluate the two contributions to ATc/T^"' in eq (15), and show how eq (17) follows from the 
criticality condition E(0,0) = fj,. 

The first integral in (15) is [/c(ek - Mc) - /c(ek)] = K^g3/2{e^'''^'') - A;r^C(|). Here the gsfzn are the Bose- 
Einstein integral functions^^, gn{x) = ^i{x^/i^), and ({n) = 5n(l). We are using the notation Jj^ = {2tt)~^ J (Pk. 
Using the Robinson expansion^^ for 5(3/2, we get 

^ [/c(ek - Mc) - /c(ek)] = ^[-2V7r(-/3cMc) + 0{M] . 
For the second integral in (17), using the second-order self-energy from eq (6), one obtains after some variable 




8 



transformations 



Jk O/Xc Jp Jq A ep + ek - eq - Cp+k-q Ag V-PcMc 

Here Cp = Cp — /Uc- In combination with the criticahty condition (17), this leads to 



The condition (17) for the transition temperature is obtained from /ic = Sr(0,0) by using the second-order pertur- 
bative result for the self-energy (e.g., obtained from eq (6) by setting k = 0): 



/(ep)/(eq 



p "^P 



Cp— q /^c 
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ln(-/3cA«c) • 



* Electronic address: masud@physics.rutgers.edu 
^ Electronic address: aiidreir@physics.rutgers.edu 

^ G. Baym, J.-P. Blaizot, M. Holzmann, F. Laloe and D. Vautherin, Eur. Phys. Journ. B 24, 107, (2001); also cond- 

mat/0107129. 

2 G. Baym, J.-P. Blaizot, M. Holzmann, F. Laloe and D. Vautherin, Phys. Rev. Lett. 83, 1703 (1999). 
^ G. Baym, J-P. Blaizot and J. Zinn- Justin, Europhys. Lett. 49, 150 (2000). 
" P. Griiter, D. Ceperley and F. Laloe, Phys. Rev. Lett. 79, 3549 (1997); 

M. Holzmann and W. Krauth, Phys. Rev. Lett. 83, 2687 (1999); 

M. Holzmann, P. Griiter and F. Laloe, Eur. Phys. Journ. B 10, 739 (1999). 
^ M. Holzmann, G. Baym, J.-P. Blaizot, and F. Laloe, Phys. Rev. Lett. 87, 120403 (2001). 
^ E. Mueller, G. Baym, and M. Holzmann, J. Phys. B, 34, 4561 (2001); also cond-mat/0105359. 

M. Bijlsma and H.T.C. Stoof, Phys. Rev. A54, 5085 (1996). 
® H.T.C. Stoof, private communication. 

^ M. Bijlsma and H.T.C. Stoof, Phys. Rev. A55, 498 (1997). 

^° A.L. Fetter & J.D. Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill, 1971. (Section 28.) 
" S. Kanno, part I in Prog. Theor. Phys. 41, 966 (1969); and part II in ibid. 41, 1145 (1969). 
12 S. Kanno, Prog. Theor. Phys. 41, 949 (1969). 

1^ M. Wilkens, F. Illuminati and M. Kramer, J. Phys. B 33, L779 (2000). 

^'^ A.M.J. Schalcel, Comment on Los Alamos preprint server, cond-mat/0004142 

15 A.M.J. Schakel, Int. Journ. Mod. Phys. B 8, 2021 (1994); 

1® "Boulevard of broken symmetries", Habilitationschrift Preie Universitat Berlin (1998); also available as preprint cond- 

mat/9805152. 
1'' H. Kleinert, cond-mat/0210162 

K. Huang, Phys. Rev. Lett. 83, 3770 (1999). 

T. Toyoda, Ann. Phys., N.Y., 141, 154 (1982). 
^° A.Z. Patashinskii and V.L. Pokrovskii, Sov. Phys. JETP 46, 994 (1964). 

21 J.D. Reppy, B.C. Crooker, B. Hebral, A.D. Corwin, J. He, and G.M. Zassenhaus, Phys. Rev. Lett. 84, 2060 (2000). 

22 M. Hasenbusch & T. Torok, Journ. Phys A, 32, 6361 (1999). 

23 H. Shi & A. Griffin, Phys. Rep. 304, 1 (1998). 

2'' For properties of the Bose-Einstein integral functions, see for example, K. Huang, Statistical Mechanics, 2nd Edition, Wiley 
1988, and R.K. Pathria, Statistical Mechanics, Pergamon 1972. Some properties (expansions) were first derived in J.E. 
Robinson, Phys. Rev. 83, 678 (1951). 



